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1. INTRODUCTION 

The speed and accuracy of maximum power point tracking (MPPT) represent important 
performances in photovoltaic (PV) systems to increase the harvested energy yield under different operating 
conditions [1]. For this purpose, many of classical and intelligent MPPT methods are recommended in the 
literature [2]—[10]. The shape and number of membership functions (MFs) used in the fuzzy MPPT affect the 
tracking performances. Meanwhile, the fuzzy MPPT method of seven triangular-shaped MFs can provide the 
best transient and steady tracking performances [11]. 

The fuzzy logic controller (FLC) can be further optimized to improve the tracking performances by 
tuning the values of MFs’ parameters using individual or combination of different intelligent techniques as 
genetic algorithm (GA), particle swarm optimization (PSO), and artificial neural networks (ANN) 
[12]-[18]. In this work, an asymmetrical FLC tracking method of seven triangular-shaped MFs is used, 
where its MFs’ values are derived based on power-voltage (P-V) characteristics under standard technical 
conditions (STC) of irradiation 1000 W/m? and cell temperature 25 °C. Moreover, the MFs of fuzzy MPPT 
method is optimized using the chimp optimization algorithm. In the optimization process, the fitness function 
is a weighted sum of a cost function under different irradiation and temperature operating conditions. Where, 
the cost function is composed of transient and steady state indices at each operating condition. The tracking 
performances of the optimized fuzzy tracking method are examined and compared with the asymmetrical and 
classical well-known perturb and observe (P&O) tracking methods in terms of rising time (t,) and energy 
yield from the PV module, and the fitness function value. 
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2. PHOTOVOLTAIC (PV) SYSTEM MODELING 

In this work, a BPSX150S PV module of 72 series solar cells is used. Under STC, this module has a 
maximum power point (MPP) of maximum power of 150 W at module voltage and current of 34.5 V and 
4.35 A, respectively [19]. Since the location of MPP is changed by changing irradiation (G) and temperature 
(T) conditions, the extracted power and energy yield from PV are subsequently changed. Consequentially, 
using a proper tracking method to maintain the MPP is an essential role to improve the PV system 
performance. Figure 1 depicts the PV system including a fuzzy MPPT controller used in this work. 

For a single-diode model, the operating point current of PV module at different module voltage (V), 
irradiation (G), and cell temperature (T) can be represented by [20], [21]: 


1 =f (V,G,T) = lpn- lo (exp (ae) = 1) = (=), (1) 


NsnKT Rsh 


where Jpn is the photocurrent or short-circuit current (Jc). Rs and Rsn are the series and parallel resistances of 
the PV cell, respectively. q is the electron charge (1.602x10-° C), K is the Boltzmann’s constant 
(1.381x10°3 J/K), n is the diode ideality factor (1.62 in this work), N, is the number of series solar cells 
constructing the PV module (72 in this work). J, is the diode’s reverse saturation current which is affected by 
the cell temperature. J, equation can be described by (2): 


I, = lor X () x exp (2 (= — »), (2) 


where T, is the temperature at STC which is 298 °K or 25 °C. E; is the band-gap energy of the semiconductor 
used in cell manufacturing. Isc which is influenced by T and G is presented as (3): 


Isc =£ (Iscr + a(T = Eo), (3) 


E 


where G; is the irradiation at STC and a is the coefficient of short circuit current temperature [20], [21]. 

Figure 2 shows the (Z-V and P-V) and the corresponding AP/AV characteristics of the PV module at 
various G and constant T of 25 °C. Where Figure 2(a) shows the Z-V and P-V characteristics while Figure 2(b) 
shows the corresponding AP/AV characteristics (absolute value). It is obvious from Figure 2(b) that the 
AP/AV is smoothly changed. Hence, it is used as a suitable input of the fuzzy MPPT controller. 
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Figure 1. PV system including fuzzy MPPT controller 
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Figure 2. At various G and constant T=25 °C: (a) I-V and P-V characteristics and (b) |dP/dV| characteristics 
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3. ASYMMETRICAL FUZZY MPPT CONTROLLER 

The Fuzzy MPPT method has promising results compared with classical tracking methods as P&O, 
InC [22], [23]. Hence, in this paper, it is used to track and maintain the MPP by continuing updating the PV 
module voltage by suitable PV voltage change (AV). Accordingly, the FLC input is the PV power-voltage 
slope (AP/AV) while the output is AV. FLC can be classified into different types according to the definition 
of setting values, shapes, and number of the used membership functions (MFs). In this paper, an 
asymmetrical FLC of seven triangular-shaped MFs is used due to its superior results compared with five bell 
and triangular MFs [11], [20]. The asymmetrical MFs of the input AP/AV is shown in Figure 3. Where 
symbols NB, NM, NS, Z represents negative big, negative medium, negative small, and zero, respectively. 
Whereas, PB, PM, PS represents positive big, positive medium, and positive small, respectively. x1, x2, x3, 
x4, x5 and x6 are setting values of asymmetrical input’s MFs. These values are determined based on the 
design strategy developed in [11]. 

Under STC, the minimum and maximum values of AP/AV are -38 W/V and 4.75 W/V, respectively 
as shown in Figure 2(b). Hence, the maximum negative value of AP/AV (x1 in Figure 3) is -38, whereas the 
maximum positive value of AP/AV (x6 in Figure 3) is 4.75. Consequently, the six setting values of MFs are 
x1=-38, x2=-25.33, x3=-12.67, x4=1.58, x5=3.16, and x6=4.75. In contrast, seven MFs of the output AV are 
shown in Figure 4, where the chosen minimum and maximum setting values are -1.5 and 1.5, respectively. 
The FLC rules can be easily concluded from observing Figure 2(b). Since the output AV is linearly depend on 
the input AP/AV. Subsequently, the seven rules depicted in Table 1 are used [11]. Furthermore, in this paper, 
an optimization method is proposed to optimize the setting values of the asymmetrical FLC MPPT method 
for improving the PV tracking performance. The next section describes the proposed tracking method. 
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Figure 3. MFs of asymmetrical FLC input AP/AV 
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Figure 4. MFs of asymmetrical FLC output AV 
Table 1. The rules of FLC 


AP/AV NB NM NS Z PS PM PB 
AV NB NM NS Z PS PM PB 


4. OPTIMIZED FLC USING CHIMP ALGORITHM 

The optimization of MFs’ setting values plays an essential role in improving the tracking 
performance. For this purpose, many algorithms can be used [15]-[18]. In this paper, the chimp algorithm is 
used to optimize the six setting values of FLC input AP/AV. In the following sub-sections, the chimp 
optimization algorithm and its details for optimizing the MFs’ setting values are presented. 


4.1. Chimp optimization algorithm 

The chimps' society is based on fission and fusion. The mixture of society in this type of society 
varies over time. Furthermore, each member of society has a unique capacity and a distinct responsibility, all 
of which may alter over time. With this in mind, this algorithm proposes the idea of separate groups, in which 
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each chimps' group attempts to learn about the search area using its unique capacity for a specific task. There 
are four types of the chimps were given the names of “driver”, “barrier”, “chaser”, and “attackers”. They 
have a variety of responsibilities in the hunting process to ensure a good hunt. The preys are pursued by the 
drivers, who do not try to catch them. Barriers erect themselves in trees to build a dam across the prey's 
escape route. Chasers run quickly after their preys in order to catch them. Finally, the attackers anticipate the 
prey's descent into the lower canopy. Attackers must be more intelligent in anticipating the prey's next 
movements. As a result, after a good search, Attackers get a bigger piece of meat as a reward. This important 
task (attacking) is linked to age, intelligence, and physical capacity. Furthermore, chimps may switch roles 
during a hunt or remain in the same role during the process [24]. To summarize, chimps' social hunting There 
are two parts to this activity: "exploration" which comprises moving, Defending and pursuing the prey, and 
"exploitation", it requires going after the prey. Figure 5 depicts the two stages. Where Figures 5(a) and 5(b) 
depict the exploration and exploitation stages, respectively. 
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(a) (b) 
Figure 5. The major stages of chimp’s operations: (a) exploration and (b) exploitation 

4.1.1. Driving and chasing the prey 

The prey is hunted during the exploration and extraction processes, as previously mentioned. 
As shown in (4) and (5) are used to model driving and chasing the prey mathematically [24]: 

Xchimp (^ + 1) = Xorey() — a.d, a = 2. fı. Rd, — a (4) 

d= |c. Xprey (N) — M. Xchimp (n)|, c = 2. Rd, (5) 
where n represents the current iteration, m is chaotic value, Xchimp (n) and Xprey (n) are the chimp and prey 
position vectors, respectively. Rdz and Rdz are the random values with the values of [0,1]. Through the 
iteration, fis nonlinearly decreased from 2.5 to 0. 
4.1.2. Exploration stage 

Two methods are designed to mathematically model chimp attacking behavior: at first, the chimps 
investigate the prey's position using driving, obstructing and pursuing, then they might be able to encircle it. 


Finally, attackers are normally in charge of the hunt. The search sometimes includes driver, barrier, and 
chaser [24]. 


dattack = |C1Xattack — Md], dbarrier = |C2Xparrier — Max, 
denaser = |C3Xcnaser — M3X|, dariver = |C4Xariver — M4 x| (6) 


X1 = Xattack — l1 lattack» X2 = Xparrier — A2dparrier 
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X3 = Xcnaser — 434 chasers X2 = Xariver — 4 ariver> (7) 


x(n) = Xz t+X24+X3+Xq 


4 (8) 


As shown in Figure 6, the chimp’s location in space of the search is modified in relation to the 
positions of other chimps. As shown, the final location of the chimp is determined at random inside a circle. 
Where a circle is identified by the positions of the attacker, barrier, chaser, and drivers. 


A- other chimps — 


Figure 6. Position update position 


4.1.3. Exploitation stage 

The chimps split up to look for prey and then they all come together to attack it. This action is 
formally modelled using the vector a, as shown in Figure 7. Where, the condition 1<|a] causes chimps to 
deviate from prey to avoiding entrapment in local optimal point and condition 1>ļa| causes chimps to 
confluence at the prey site to get global optimal point. 
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Figure 7. The parameter “a” effecting on the updating operation 


4.1.4. Improved the exploitation stage 

As previously mentioned, the social reward in chimp society is based on hunting meat. In the final 
phase, chimps are forced to give up their hunt behavior by attempting to steal the meat. Own hunting 
responsibilities as a result, they attempt in a haphazard manner to steal hunting meat for social purposes. 
Chaotic maps can be used to model the chimps' chaotic behavior in the final stage. The updating operation 
can be modeled as [24]: 
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Xprey(n)-axd ifu<0.5 


Chaotic value if u 2 0.5 ©) 


Xchimp (n+1)= { 


where u is a value between zero and one. 


4.2. Optimization of MFs’ setting values using chimp algorithm 

The chimp algorithm is used to optimize the MFs’ setting values of FLC input AP/AV: x1, x2, x3, 
x4, x5, and x6 by optimizing the intervals between setting values: xxl, xx2, xx3, xx4, xx5, and xx6. 
Meanwhile, the optimization process can be defined by maximizing the fitness function f (xxl, 
xx2, XX3, xx4, xx5, and xx6). Where x(1)=-xx(1)-xx(2)-xx(3), x(2)=-xx(1)-xx(2), x(3)=-xx(1), x(4)=xx(4), 
X(5)=xx(4)+xx(5), x(6)=xx(4)+xx(5)+xx(6), based on the constraints: 


O<xx 1+xx2+xx3<NEG max, 

0<xx4+xx5+xx6 <POS max, 
where the selected values of NEG max and POS max in this work are 50 and 10, respectively. 

The goal of optimization is maximizing the energy yield from the PV module and minimization the 
rising time at different irradiation and temperature operating conditions. Subsequently, in this paper, the 
chosen fitness function “Fitness” is a weighted sum of a cost function at five different operating conditions of 
G and T (in Iraq). Hence, the fitness function can be defined as (10): 


Fitness = X}; Cost(i) w(i) (10) 


where Cost(i) is the cost function which is composed of transient and steady-state response indices at the 
operating condition ¿i as (11): 


tf 
7 Íe- V(t)-I(t)dt 


Cost(i) = fos (=) +0. | x 100% (11) 
f 


tf 
Se Pmppat 


where ¢ is final time of simulation [11]. Whereas, w(i) is the weight of the operating condition i. The best 
fitness function value is 1. 

In this work, the chosen values of weights are 0.25, 0.166, 0.084, 0.33, and 0.17 corresponding to 
the five chosen operating conditions (in Iraq) defined by G(i) of 400 W/m’, 600 W/m?, 1000 W/m?, 
800 W/m?, 1000 W/m”, and T(i) of 25 °C, 35 °C, 25 °C, 45 °C, and 45 °C, respectively. According to each set 
of setting values, the fuzzy MPPT is simulated at five different operating conditions to calculate the 
corresponding fitness function value. Finally, the best solution is chosen. 


5. RESULTS AND DISCUSSION 

During the optimization and test phases, the chosen initial operating PV voltage is 1 V (at left of 
MPP). The performance of the optimized tracking method is examined by MATLAB simulation for 30 s at 
different operating conditions. Furthermore, results of the proposed method are compared with the classical 
P&O method and asymmetrical FLC method (explained in section 3) in terms of rising time (f,) and energy 
yield from the PV module. Where the energy yield in (Wh) can be expressed by [7]. 


tf 
So VOMoat 
3600 


Energy yield = (12) 

According to the optimized FLC by chimp algorithm, the optimized solution is xx1=5.3364, 
xx2=4.8655, xx3=6.4194, xx4=0.0048, xx5=0.9542, and xx6=0.9398. Consequently, the MFs’ setting values 
of AP/AV are x1=-16.6213, x2=-10.2019, x3=-5.3364, x4=0.0048, x5=0.959, and x6=1.8988. Hence, 
Figure 8 is shown the corresponding MFs of asymmetrical and optimized FLC tracking methods are shown in 
Figures 8(a) and 8(b), respectively. 

Figure 9 is shown at STC and by using the different tracking methods, Figures 9(a) and 9(b) show 
the PV power and energy yield, respectively. Where maximum output power and energy can be harvested 
from the used PV module under this condition are 149.9873 W and 1.2499 Wh, respectively. In the same 
context, Table 2 illustrates the comparative results of the tracking methods under STC. In contrast, Figure 10 
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shows the power as shown in Figure 10(a) and energy yield in Figure 10(b) curves, respectively at operating 
condition of G=400 W/m? and T=25 °C. Where the maximum output power and energy can be harvested 
from the module under this condition are 56.8948 W and 0.4741 Wh, respectively. Whereas, the comparative 
results of the tracking methods under this operating condition are illustrated in Table 3. Since the MFs’ 
setting values of asymmetrical FLC method is derived under STC, hence its tracking performances are 
obviously close to the tracking results of optimized method at STC as shown in Figure 9 and Table 2. 
Whereas the asymmetrical FLC method failed to satisfy the optimal results under the other operating 
condition of G=400 W/m? and T=25 °C as shown in Figure 10 and Table 3. 
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Figure 8. MFs of AP/AV: (a) asymmetrical FLC MPPT and (b) optimized FLC MPPT 


In contrast, the optimized tracking method can successfully track the MPP with less rising time 
compared with asymmetrical method. Consequently, a more energy is harvested under different conditions as 
shown in Figures 9 and 10. Where in the case of G=400 W/m? and T=25 °C, energy of 0.4581 Wh with rising 
time of 2 s are harvested by using the optimized method, in contrast, energy of 0.4386 Wh, 0.4515 Wh, and 
0.4299 Wh with rising times of 4.4 s, 0.8 s, 5.3 s are harvested using asymmetrical, P&O (at AV=3.5 V), and 
P&O (at AV=0.5 V) tracking methods, respectively. 


150 7 ag 1 a P&O (AV=0.5 V) 
a) i a int na ny tity Z 4.25 ~~ ~ PO (AV=3.5 V) 
z R E+ 
3 ini ORHAN = == Asymmetical FLC Z 
is oO Optimized FLC 
© > 1 A 
= 100 x — — —E ideal 
[0] 
a S 0.75 
E i 
3 P&O (AV=0.5 V) 2 o5 
= 50 — — — P80 (AV=3.5 V) = 
A Asymmetrical FLC g 
Optimized FLC > 
oO 


— — = P ideal 


0 5 10 15 20 25 30 0 5 10 15 20 25 30 


Figure 9. PV power and energy yield under STC: (a) power and (b) energy yield 


Moreover, the performance of the optimized method is compared with the other tracking methods in 
terms of fitness function value. Table 4 depicts the comparison results. In future the system can be improved 
using different objective function as in [25]-[27]. 
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Table 2. Comparative results of tracking methods under STC 


Tracking Method Rising Time (s) Energy Yield (Wh) 
P&O (at AV=0.5 V) 5.6 1.1269 
P&O (at AV=3.5 V) 0.8 1.1525 
Asymmetrical FLC 2.2 1.2051 
Optimized FLC by chimp algorithm 2.1 1.2055 
60 0.5 r r r 
z P&O (AV=0.5 V) 4 
eal = — — —P&0 (AV=3.5 V) 
= a 04 |= = Asymmetrical FLC 
5 io g Optimized FLC 
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Figure 10. PV power and energy yield at G=400 W/m? and T= 25 °C: (a) power and (b) energy yield 


Table 3. Comparative results of tracking methods at G=400 W/m? and T=25 °C 


Tracking Method Rising Time (s) Energy Yield (Wh) 
P&O (at AV=0.5 V) 5.3 0.4299 
P&O (at AV=3.5 V) 0.8 0.4515 
Asymmetrical FLC 4.4 0.4386 
Optimized FLC by chimp 2 0.4581 
algorithm 


Table 4. Fitness function values of tracking methods 


Tracking Method Fitness 
P&O (at AV=0.5 V) 0.9468 
P&O (at AV=3.5 V) 0.9558 
Asymmetrical FLC 0.9704 


Optimized FLC by chimp algorithm 0.9802 


6. CONCLUSION 

In this paper, an asymmetrical fuzzy MPPT method of seven triangular MFs is used. Where, its 
setting values are derived under STC. Furthermore, the fuzzy MPPT method is optimized by using the chimp 
algorithm taken into consideration different operating conditions included in the chosen fitness function. It is 
demonstrated from the simulation results in section 5, the asymmetrical FLC method cannot satisfy the 
optimal tracking performances at all test operating conditions. Nevertheless, its tracking effectiveness is 
conspicuous at STC. In contrast, the optimized tracking method by chimp algorithm can successfully 
investigate the optimal tracking performances at all test conditions. Moreover, it is evident that the optimized 
method has a largest fitness of 0.9802 compared with asymmetrical, P&O (at AV=3.5 V), and P&O (at 
AV=0.5 V) methods of 0.9704, 0.9558, and 0.9468 fitness values, respectively. Hence, it can be concluded 
that the optimized tracking method can successfully enhance the transient and steady state tracking 
performances for the PV systems. As future work, a fuzzy type-2 and another intelligent hybrid MPPT 
algorithms have to be explored and examined to conclude the best tracking algorithm for improving the 
overall PV performances. 
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